The role of lubricant molecular shape in microscopic friction 
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With the help of a simple two-dimensional model we simulate the tribological properties of a thin 
lubricant film consisting of linear (chain) molecules in the ordinary soft-lubricant regime. We find 
that friction generally increases with chain length, in agreement with their larger bulk viscosity. 
When comparing the tribological properties of molecules which stick bodily to the substrates with 
others carrying a single sticking termination, we find that the latter generally yield a larger friction 
than the former. 

PACS numbers: 81.40.Pq; 46.55.-fd; 62.20.Qp 



I. INTRODUCTION 

Understanding sliding friction between substrates sep- 
arated by a thin lubricant film is technologically crucial 
but also rich physically [l|, 0]. Without lubricant, it is 
known that the lowest friction can be achieved for a con- 
tact of two hard crystalline solids, when the atomically 
flat surfaces are incommensurate (e.g., see [2] and ref- 
erences therein). When the surfaces are separated by a 
lubricant film, the lowest friction is achieved if the film is 
solid and crystalline but incommensurate with the sub- 
strate surface. Such a situation appears in the so-called 
"hard" lubricant system [^l , where the intermolecular in- 
teraction strength within the lubricant, Vh, is stronger 
than the lubricant-substrate interaction, Vsi. In that 
case, the shape of lubricant molecules is not too relevant, 
because friction is determined uniquely by the substrate- 
lubricant interface structure. 

Conventional lubricants, however, belong typically to 
the opposite "soft" lubricant type, with Va < Vsi. Here 
each substrate is covered by a glued lubricant monolayer, 
protecting it from wear. Sliding occurs deeper inside the 
lubricant film, where below the film melting point a typ- 
ical melting- freezing mechanism operates: the lubricant 
locally melts at slip (i.e., during sliding) and re-solidifies 
at stick Q. For a thick lubricant film above its melting 
point, e.g., and for thicknesses larger than say ten molec- 
ular diameters, friction is proportional to the bulk viscos- 
ity of the lubricant. Therefore, thick fluid lubricants con- 
sisting of simple spherical molecules or individual atoms, 
characterized by a low bulk viscosity, would provide the 
best frictional properties. Indeed, experiments show that 
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the bulk viscosity for linear molecules grows with the 
molecule length as was shown in Ref. [J|. One could 
then expect that the friction coefficient of chain lubri- 
cant molecules should grow with their molecular length. 

However, it is well known that long molecules often 
make better lubricants. The reason for this is the squee- 
zout effect: at high load, when the regime of boundary lu- 
brication is approached and the film is locally extremely 
thin, at stick the lubricant may be squeezed out in cor- 
respondence of the substrate asperities. Direct contact 
of unlubricated asperities leads to a high static friction 
force fs and generally to substrate wear. Sivebaek et 
al. 4] showed that with the increase of the alkane chain 
length, transitions from Ni to TV; — 1 lubricant layers oc- 
curs at higher pressure. Longer lubricant molecules are 
harder to squeeze out, and this is why they provide bet- 
ter lubrication properties. Following that line, Ref. [3| 
studied the "totally-glued" case, where all atoms of the 
lubricant molecule interact strongly with the substrates. 
Longer molecules adhere more to the substrates, making 
it harder for the squeeze-out of the last lubricant layer 
to nucleate. In agreement with that, longer alkanes per- 
form better as boundary lubricants than shorter ones, 
because they prevent more effectively the appearance of 
cold-welded junctions and the resulting wear. 

Tribological properties of lubricants constructed of lin- 
ear (chain) molecules of length L, again for the totally- 
glued system, were also studied with the help of molec- 
ular dynamics (MD) simulation by He and Robbins [H]- 
They found that at the small coverage 9 = M/Ng = 1/8 
(the number of lubricant molecules M divided by the 
number of substrate surface atoms Ng) the dependence 
of the static friction on chain length (for L — 1,3,6) is 
small, while for a higher coverage, M/Ng — 1/2, chains 
with L = 3 and 6 show nearly the same behavior, while 
"monoatomic molecules" (L = 1) display a dramatically 
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FIG. 1: (Color online) A schematic cartoon of our simula- 
tion model. This snapshot refers to a configuration where the 
Ns — 31-atoms substrates (blue) squeeze M = 35 head-glued 
molecules, each made of L = 5 atoms (red). 



reduced friction (about one fourth). 

The aim of the present work is to study the role of 
molecular length in detail for a wide range of model pa- 
rameters. Moreover, we wish to compare a totally-glued 
lubricant (with all monomers of the lubricant chain stick- 
ing to the substrates) to a "head-glued", with a single 
(head) atom of the molecule strongly attached to the 
substrate, while the other atoms (molecular tail) interact 
more weakly with the substrates. Our aim is to clarify the 
role of the length of the lubricant molecules in friction, 
in particular to find what length of lubricant molecules 
could provide best tribological characteristics for a thin 
(few molecular layers) soft lubricant film. We will use 
the following criteria to define good tribological proper- 
ties: (i) a better lubricant provides lower values for both 
the static friction force fs and the kinetic friction force 
fk, and (ii) a better lubricant is characterized by lower 
values for the critical velocity Vc of the transition from 
stick-slip to smooth sliding (or at least to provide a less 
irregular motion in the stick-slip regime). 

With the aim of addressing general trends, we choose 
to explore a simplified two-dimensional (2D) minimal 
model which in our view is sufficient to catch the physics 
of the problem, with qualitatively correct trends. In re- 
turn for the model's simplicity, we will be able to span 
a large choice of structures and parameters. On top of 
that, the 2D geometry provides an easier visualization of 
all processes inside the lubricant. The 2D model does not 
describe realistically the lubricant squeezout [5], which 
requires a full three-dimensional (3D) modeling, such as 
that employed in MD simulations of Ref. [4] . Therefore 
we will not address squeezout properties any further. 

The paper is organized as follows: Sec. |TT] introduces 
our minimal model; Sec. IIIII provides the mathematical 
and numerical details of the model's solution; Sec. IIVI 
contains a detailed analysis of the results, which are fi- 
nally discussed in Sec. |Vl 



II. THE MODEL 

We use a 2D model, where point particles, representing 
individual atoms or monomer units of larger molecules, 
can move in two dimensions x and z, where x is the slid- 
ing direction and z is perpendicular to the substrates. 
The substrates are modeled as two rigid chains of Ng 
atoms and lattice spacing Rg, so that the total mass of 
each substrate is Ngiris and the system size in the sliding 
direction is = NsRg. We apply periodic boundary 
condition along the x direction. The bottom substrate 
is fixed at x = z = 0, while the top one is free to move 
in both X and z directions. The top substrate is pressed 
toward the bottom one by a (— z)-directed constant load 
force Fi = Ngfi, and is driven in the x direction at ve- 
locity Vs through an attached spring of elastic constant 
Ngks, as sketched in Fig. [TJ The spring force F, which 
works against the friction force, is monitored during sim- 
ulation (throughout the paper we report the force per 
substrate atom / ~ F/Ng). Thus, our model is a 2D 
variant of a typical experimental setup used in tribol- 

ogy[l|,i. 

Between the substrates we insert N lubricant atoms 
of mass mi. How eventually these N atoms are to be 
lumped together to form linear molecules will be detailed 
below. Independently of that, all atoms interact via pair- 
wise 12-6 Lennard- Jones (LJ) potential 



14j(r) = T4 



Rao.' 



(1) 



Here a, a' = s or I for the substrate or lubricant atoms 
respectively. Thus, the lubricant-lubricant interaction is 
characterized by the parameters Vu and Rh, while the 
lubricant-substrate interaction by Vsi and R^i (direct in- 
teraction between the top and bottom substrates is omit- 
ted). Throughout the paper we will use dimensionless 
units, where — mi = 1, Ru — 1, and the energy pa- 
rameters Vaa' takes values around Vaa' ~ 1- The relative 
strength of the energy parameters Vu and Vgi determines 
the low-temperature lubricant behavior The lubri- 
cant is "hard" (i.e., at low temperature it remains solid 
at slip) when Vu > Vsi, and "soft" (i.e., the film melts 
at least locally during sHps) for Vu <C Vsi- Because a 
2D model cannot reproduce even qualitatively the true 
phonon spectrum of a 3D system, and because moreover 
frictional kinetics is generally diffusional rather than in- 
crtial, we use Langevin equations of motion with a Gaus- 
sian random force corresponding to a given temperature 
T, and a damping force 



fr,,x = -miri{z)x - miri{Z - z){x - X) . 
fn,z = -~miri{z) z - mir]{Z - z) {z - Z) . 



(2) 



Here the coordinates of a generic lubricant atom 

and X, Z is the center-of-mass coordinate of the top sub- 
strate. The viscous damping coefficient is assumed to 
decrease exponentially with the distance from the corre- 
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spending substrate, 



ri{z) = 770 [1 - tanh(z/z<j)] 



(3) 



where we typically use 770 = 1 and — 1. This form 
is meant to mimic the Joule heat dissipation, which can 
only take place through the substrates, and therefore is 
mostly effective at the interfaces. The Gaussian width of 
the Langevin random force equals [2ri{z)miTY^'^ 

We compare the main tribological properties of lubri- 
cants consisting of molecules of different lengths, from in- 
dividual (atom-like) monomers to linear chain molecules 
composed of L atomic/mo nomeric units, so that the 
number of lubricant molecules is M = N/L. To build 
chain molecules, we follow the procedure due to Rob- 
bins et al. [1, Q. The lubricant is described with the 
Kremer-Grest bead-spring model fej, which yields realis- 
tic dynamics for polymer melts [lO| . While all monomers 
within a molecule still interact with each other via the 
LJ potential, Eq. ([1]), adjacent monomers inside a given 
molecule (chain) interact via an additional ( "FENE" ) po- 
tential [ii 



for r < i?cH 
for r > i?cH, 



(4) 



where i?cH = 1.336 i?;; and Vch = 33.75 V/; (these pa- 
rameters are the same as in Ref. 3] ) . This additional po- 
tential has the main effect of preventing molecular break- 
ing. Besides that, it is fairly small at the equilibrium 
distance r ~ Ru of Vlj, and causes a bond-length con- 
traction to -Rmol — 0.856i?;i. 



III. SIMULATION TECHNIQUE 

Our frictional simulations use periodic boundary con- 
ditions throughout, a choice which as usual minimizes 
size effects. Results are generally dependent on the total 
number of lubricant atoms N. If N does not match ex- 
actly the exact number of atoms required for an integer 
number of closely packed layers, then the extra inter- 
stitial atoms or vacancies constitute structural defects. 
In real life the presence of these defects is of course the 
rule rather than the exception, but in the large interface 
size of an actual experiment their importance is more 
marginal than in our finite-size simulation. To reduce er- 
rors due to these misfit defects, we use a geometry with a 
slightly curved top substrate, whose z-coordinate varies 
along the x direction by 
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as illustrated in Fig. [T] 

In the simulations to be presented below, each sub- 
strate contains Ns — 62 rigid atoms and the substrate 
lattice constant is chosen a.s Kg = 2/3, moderately 



"incommensurate" with the lubricant equilibrium inter- 
atomic distance Ru = 1, a feature not atypical for lubri- 
cant/substrate interfaces. 

Accordingly, the system size is = 41.3. This fits up 
to < 50 atoms (in the atomic lubricant case) in one 
monolayer film under a standard applied load fi — I. 
The N = 100 atoms therefore complete two monolayers 
(one glued to the bottom substrate, the other to the top 
substrate). Atoms in excess of 100 pile up in between, 
giving rise to a third, then fourth, etc., layer, consisting 
of up to 50 atoms, for the typical load conditions of the 
simulations. 

To guarantee an efhcient melting-freezing of the soft 
lubricant in correspondence to the stick-slip motion for 
the "incommensurate" lubricant-substrate interface, we 
require a rather strong inequality Vsi ^ Vu . Specifically, 
we choose Vu = 1/9 and Vsi — 5. In detail, we set the 
lubricant-substrate interaction energies to be Vsi — 5 in 
the "totally-glued" case, while for the "head-glued" sys- 
tem, the large interaction Vsi = 5 applies exclusively to 
the "head" molecular atom, while other atoms in the rest 
of the molecule (the "tail") interact with the substrates 
with the much weaker V^i = 1/9. With these parameters, 
the melting temperature of the lubricant is of the order 
of T « 0.1; it is a generally decreasing function of the 
lubricant layer thickness, but an increasing function of 
the chain length and applied load [ll|. 

The load force used is typically fi — I per substrate 
atom. Assuming an energy scale 9Vu — 0.1 eV and inter- 
atomic distances Ru — 0.3 nm (typical of Van der Waals 
intermolecular interactions) , in a hypothetical 3D geom- 
etry such load would correspond to a pressure of almost 
1 GPa. 

The top substrate is driven at an average velocity 
Vs through a relatively stiff spring of elastic constant 
ks — 0.1. For the remaining parameters, we take 
Rsi = Rs and z^ = Rs- We run all simulations start- 
ing from well-relaxed and annealed initial configurations, 
and exclude from averaging any initial transient leading 
to the steady dynamical state. In more detail, we solve 
the Langevin equations of motion using a fourth-order 
Runge-Kutta method with a time step At = ro/128, 
where tq = 2i:Ru{72 ■ 9Vu/mi)^^^'^ is an average char- 
acteristic period of the LJ interatomic potentials. The 
initial configuration is realized by a careful annealing cy- 
cle up to the high temperature T ~ 2, followed by a 
"relaxation" run, typically of 4096 tq or longer, at the 
appopriate Langevin temperature T and sliding velocity 
Vs , before measuring the friction force by averaging over 
a run typically of the same duration. The final configu- 
ration of every run is saved: this allows us to restart the 
run when we need to improve accuracy. To estimate error 
bars of various quantities, we split the measuring trajec- 
tory into 30 pieces of equal duration, and then estimate 
the variance of the global average through the standard 
deviation of the averages carried out over each individ- 
ual piece. To verify the results, for some parameters we 
also made runs for total integration times and system 
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FIG. 2: (Color online) Comparison of the stick-slip motion 
[vs = 0.03, panels (a,c,e)] and the smooth sliding [vs = 0.1, 
panels (b,d,f)] for M = 40 totally-glued three-atomic {L — 3) 
lubricant molecules at T = 0. Panels (a) and (b): typical 
time dependence of the friction force f{t); panels (c) and (d): 
histogram P{f) of the friction force f{t); panels (e) and (f): 
depth distribution P(C)- The top substrate is curved, see 
Fig. m and the parameters are the following: Vu — 1/9 and 
Ru = 1, Vsi = 5 and R,i = R^, fi = 1, fc, = 0.1, Rs = 2/3, 
and Ns = 62. 



sizes 2-^4 times larger to check that the results do not 
change. 

As in typical tribological systems, we observe a tran- 
sition from stick-slip to smooth shding with increasing 
driving velocity. An example is shown in Fig. [21 In the 
smooth sliding regime, the kinetic friction force can eas- 
ily be measured as the average of the spring force f{t). 
To evaluate the static friction force fs in the stick-slip 
regime, one can calculate the histogram P(/) of f{t) and 
then select fs as the point where P{f) decays below a 
given level, e.g., P{fs) — 0.5 maxP(/) (see Fig. 
A harder problem is that of determining the precise 
transition from stick-slip to smooth sliding and its crit- 
ical velocity Vc, since this transition is often a smooth 
crossover. For this purpose, we calculate the distribution 
(histogram) P(C) of an auxiliary variable Ci defined as 
C(i) — [df{t)/dt] I (ksVs). In smooth sliding the function 
P{C) decays monotonically from a maximum at C ~ 
regime (see Fig. [2f). In contrast, for stick-slip, P{C) 
should show a maximum at some < < 1 (Fig. 
Indeed, during sticks the force grows as f{t) ~ kgVst [so 
that P{Q = 6{( — 1)], and it drops quickly to substan- 
tially smaller values during short slip events. Thus, this 
method allows us to evaluate Vc fairly accurately. 



IV. RESULTS: TOTALLY-GLUED VERSUS 
HEAD-GLUED MOLECULES 



For the chosen set of parameters, the transition from 
stick-slip to smooth sliding occurs at a velocity Vs < 1] 
therefore, we select = 1 to represent the smooth- 
sliding regime, where we calculate the kinetic friction 
force. Figure [3] summarizes the simulation results for 
T = 0, reporting the time-averaged kinetic friction force 
fk as a function of the total number of lubricant atoms, 
which is basically proportional to the lubricant layer 
thickness, for different molecular lengths L. 

The monoatomic lubricant (diamonds in Fig. [S]) forms 
two complete monolayers for A^ = 100 and three mono- 
layers for N — 150. For N > 150, one ordered layer sticks 
to the bottom substrate and another ordered monolayer 
to the top substrate, while the layer(s) in between are 
liquid during sliding. For the monatomic lubricant, the 
static friction force (not shown) is larger than the ki- 
netic friction by a factor of 2 to 4; the transition from 
stick-slip to smooth sliding occurs at Vc ^ 0.03. Fric- 
tion generally decreases as the number of lubricant atoms 
N - and thus the lubricant thickness - grows, but the 
dependence fkiN) cannot be described by the simple 
laws of viscous-friction flow. At lower lubricant thick- 
ness 100 < N < 150, when the middle layer is incom- 
plete, friction is slightly lower than for three complete 
monolayers, see Fig. [31 When the lubricant amount de- 
creases further, 70 < A^ < 100, there remain only two, 
generally incomplete, lubricant layers, each glued to the 
nearest substrate (e.g., A^ = 100 atoms form exactly two 
complete lubricant layers, which slide slowly, thanks to 
incommensurability, over the substrates to which they 
are attached) . At even smaller values of A^ we have one 
monolayer or eventually just a monolayer island in the 
narrow gap between substrates, and at this point both 
the friction and the critical velocity Vc grow rapidly when 
A^ decreases. 

For all chain molecules (polymer lubricant) L > I, fric- 
tion is always larger than for the monomer lubricant. The 
kinetic friction force increases with the chain length L, 
if the total number of atoms N is kept fixed. In gen- 
eral, regardless of L, friction decreases roughly mono- 
tonically with increasing TV, the lubricant thickness. In 
Fig. m we re-organize the kinetic friction fk data of Fig. [31 
but as a function of the number of molecules M (only 
points corresponding to three or more molecular layers 
are shown). In the standard viscous flow of a fluid with 
a viscosity rj^ the friction force fk should depend on the 
thickness d = (Z) of the liquid layer {d oc N oc M) as 
fk — TjvsR^/d. Unlike expectations for a fluid lubricant 
at high temperature, at T = the log-log plot of the 
function fk{M) suggests that the simulation data may 
be (crudely) flt by power laws fk ^ 0.18 (M/ATJ-^/^ ^j. 
the totally-glued case, and fk « 0.26 {M/Nsy^'^ for the 
head-glued system. Accordingly, the function fk{M/Ns) 
displays a nontrivial apparent "universal scaling" with 
thickness, at least within the nanoscale thickness range 
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FIG. 3: (Color online) Smooth sliding (T = Q, Vs = 1) time-averaged kinetic friction force fk as a function of dimensionless 
coverage N/Ns (bottom scale), as well as total number of atoms A'^ = ML (top scale) for linear lubricant molecules of different 
lengths L = 1 (black diamonds), 3 (red circles), 5 (magenta up triangles), and 7 (blue down triangles), for either the totally- 
(open symbols and dotted lines) or head-glued (solid symbols and curves) cases. The top substrate is curved according to 
Eq. ((51 (see Fig. [T]), and the other parameters are the following: Vu — 1/9 and Ru = 1, Vsi = 5 or 1/9 and Rsi = Rs, fi = 1, 
ks = 0.1, Rs = 2/3, and iV^ = 62. 



considered here. The slower decay of friction compared 
to may be attributed to a non-uniformity of local 
heating due to sliding- induced local shear. At fixed driv- 
ing velocity, the local shear decreases as thickness in- 
creases, whence the local viscosity involved increases, and 
so does frictional dissipation, causing a positive devia- 
tion from the fluid-like behavior. In other words, 
the higher average kinetic energy, and therefore higher ef- 
fective temperature, produces a smaller viscosity rj p^ . 
thus a smaller friction for a thinner lubricant layer than 
for a thicker one, moderating the power-law decay of fk 
in Fig. 3) For a much thicker lubricant, a T — simula- 
tion would eventually show lubricant solidification with 
slip movements associated to local melting, even at the 
very high sliding velocity considered here. 

Figure reports similar calculations carried out at a 
finite temperature T = 0.2, compared to which the ex- 



tra kinetic energy induced by shearing at sliding velocity 
Vs = 1 becomes neglible when the lubricant thickness 
exceeds about 100 i?/;. When the lubricant thickness ex- 
ceeds such values we do observe a crossover to the ex- 
pected scaling fk oc {M/Ns)~^ oc , as soon as the 
shearing kinetic energy becomes smaller than the ther- 
mal fluctuations imposed by the finite T = 0.2. Note 
that small finite-size effects induced by the periodic simu- 
lation box tend to delay slightly the onset of the fk oc 
constant-viscosity regime. 

Analyzing the simulation movies for the totally-glued 
lubricant with L > 1 and large enough N we observe two 
complete layers sticking to the corresponding substrates 
and a liquid lubricant between these layers. However, 
not all boundary molecules are glued completely, some 
of them are glued partially, i.e. only a part of a molecule 
sticks to the substrate while another part is detached 
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FIG. 4: (Color online) Smooth sliding. Kinetic friction force 
fk as a function of the number of lubricant molecules M (of 
different lengths L) in a log-log scale. Parameters and nota- 
tion are as in Fig. [3] Dashed lines describe power-law fits. 
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FIG. 5: (Color online) Kinetic friction force fk as a function of 
the lubricant thickness d = (Z) (roughly proportional to A*') 
in log- log scale. The parameters and notation are as in Fig.O 
except for T = 0.2. The dotted line sketches the fk <x 
thick-lubricant regime characterized by a constant viscosity 
rj. The dashed line with error bars illustrates the effect of the 
finite simulation size on fk in the L — 1 case. 



and points into the lubricant bulk. It is precisely these 
"free tails" that make the friction higher as compared 
with the monoatomic case. At lower lubricant thickness 
(M < 50 for the L = 3 case, M < 35 for the L = 5 
case, and M ^ 25 for the L = 7 system) the middle 
layer is incomplete: its molecules group together at the 
narrowest intersubstrate gap and just ahead of it (jam- 
ming/plowing effect). With a further decrease of lubri- 
cant thickness, sliding takes place at a single interface 
only, namely between the lubricant layers which are at- 
tached to bottom and top substrates. In this case, we 
observe a slow "caterpillar" motion of glued molecules 
relative the substrate, but this effect is weak because of 
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FIG. 6: (Color online) z-distribution of the average den- 
sity (circles) and atomic velocities for (a) the totally-glued 
molecules (open symbols) and (b) the head-glued (solid sym- 
bols) system with L — 5 and M = 70 (A^ = 350; other param- 
eters as in Fig. |3}. Triangles stand for the total velocity 
diamonds for the a;-component of the velocity. In the right- 
side panels z is measured with respect to the (fluctuating) top 
layer position Z, and then translated by its time average (Z). 
The apparent discontinuity at the center is therefore not real. 



the strong lubricant-substrate interaction. At even lower 
lubricant thickness occasional molecules may be glued 
to both surfaces simultaneously, and in such cases the 
molecular end sticking with fewer atoms slides along the 
corresponding substrate. 

We come now to analyze the head- glued lubricant, com- 
paring its behavior to the totally-glued one described 
above, all other parameters being the same. Firstly, the 
head-glued system (solid symbols in Figs.[2]and|3]) display 
a systematically higher friction than the totally-glued lu- 
bricant (open symbols). As in the totally glued case, the 
surfaces are covered by lubricant layers. Again, the lay- 
ers are orientationally ordered: in this case the molecular 
heads are attached to the corresponding surface, while 
the tails are directed outside the surface, and inclined in 
the sliding direction, like combed hair. Also, the sliding 
interface is more "irregular" , than in the totally-glued 
system. Therefore, one may expect a smoother frictional 
motion even in the stick-slip regime at low driving veloc- 
ities. It is important to note that contrary to the fully 
glued case, in the head-glued system we find essentially 
no jamming effects even at low lubricant thickness. 

Figure [S] depicts the different atomic density p{z) and 
velocity distributions Vx{z) and \v{z)\ across the lubri- 
cant for the totally- and head-glued molecular lubricants. 
In both cases, one layer sticks to the substrate, but the 
thickness of this layer is ~ Ru in the totally-glued system 
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FIG. 7; (Color online) z-distribution of the total energy losses 
rj{z) p{z) (z) (up triangles) and the x- losses 77(a) p(«) (z) 
(down triangles) for the totally-glued (blue open triangles) 
and head-glued (red solid triangles) lubricants for the model 
parameters as in Fig. [S] The z-integrated energy loss E is 
indicated for each curve. 



and much thicker, ^ L Rnwi, similar to the molecular gy- 
ration radius, in the head-glued case. Accordingly, while 
the sliding profile Vx{z) is fairly linear in the totally-glued 
case, in the head-glued case the shearing dissipation con- 
centrates mostly in a bulk lubricant region remote from 
the substrates. 

To elucidate the different tribological properties of the 
two kinds of molecular lubricants, in Fig. [7| we analyze 
the distribution of energy losses (dissipated power) due to 
sliding, with the method described in Ref. [3| ■ The losses 
occur mainly within the first (glued) lubricant layer, as a 
consequence of the model Eqs. ^ and ([3]), mimicking the 
capability of different lubricant layers to transfer energy 
into phononic/electronic degrees of freedom of the sub- 
strates. Importantly, we find that the main losses come 
from atomic vibrations in the transverse (z) direction, by 
far dominant over those coming from the motion along 
the driving (x) direction. This observation explains why, 
even though head-glued molecules displace sliding away 
from the substrates surfaces where most dissipation takes 
place, the head-glued system leads to a higher friction: 
transverse vibrations of larger amplitude penetrate bet- 
ter near the substrates (where most dissipation occurs) 
through head-glued molecules than through totally-glued 
ones. 

Our results suggest that a totally-glued lubricant is 
more effective than a head-glued one, because it provides 
lower kinetic friction. According to Sivebaek et al. [1], 
totally-glued molecules with their stronger sticking to the 
substrate, should also be more stable against squeezout 
of the contact area at high load. Unfortunately, as was 
said in the beginning, in our simplified 2D model we can- 
not study quantitatively the squeezout process. Squee- 
zout involves nucleation of a small [N — 1) negative is- 
land or "crater" and its subsequent expansion, a process 
which also involves elasticity of the substrates, and which 
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FIG. 8: (Color online) Dependence of the critical transition 
velocity Vc (circles) between stick-slip and smooth sliding, 
of the static friction force (red down triangles) and kinetic 
friction force (blue up triangles) at the crossover velocity 
Vs = Vc upon the dimensionless coverage 6m (bottom scale), 
or the number of lubricant molecules M (top scale) for the 
L = 5 chain molecules, for (a) totally-glued and (b) head- 
glued molecules. Simulation parameters same as in Fig. [S] 



is very characteristically three-dimensional [6]. We can, 
however, analyze in detail the dependence of frictional 
forces on the number of lubricant molecules M. Results 
for the L = 5 lubricant molecules are presented in Fig. [51 
It is useful to define the maximum number of molecules 
M2 which can arrange themselves exactly to form the 
two interface molecular layers (one attached to the bot- 
tom substrate and the other to the top substrate), so 
that sliding takes place in the remaining lubricant just in 
between these two layers. For the totally-glued lubricant 
M2 w 2Lx/[Rii + {L - l)i?inoi], and for the head-glued 
lubricant M2 < 2Lx/Ru- M2 depends on L and on the 
load force, so that it takes an analysis of the MD trajec- 
tories to define its precise values. In the current model 
(with = 41.3, Ns = 62, L = 5) we find M2 = 20 for 
the totally-glued system and M2 = 58 for the head-glued 
case. Based on the definition of M2 we can introduce the 
dimensionless coverage 9m = 2M/M2, so that 9m = 2 
corresponds to precisely two molecular layers. Figure [8] 
shows that, as expected, both the friction force and the 
critical velocity Vc increase when the coverage decreases 
below the two-layer value. The rise of friction below AI2 
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FIG. 9: (Color online) Dependence of the static (red down 
triangles) and kinetic (blue up triangles) friction forces at 
Vs ~ Vc and of the critical velocity (circles) on the length 
L of lubricant chains for the coverage Qm of two molecular 
layers, in the cases of (a) the totally-glued and (b) the head- 
glued system. Other parameters same as in Fig. |31 



= 2 is fairly slow in the head-glued system, but is very 
sharp in the totally-glued case. As is well known [1], in 
real 3D system the lubricant is squeezed out from the 
contact areas at the asperities, where about 2 — 3 lubri- 
cant boundary layers are left, harder to remove. When 
however the squeezing force is pushed further leaving less 
than two layers, the head-glued lubricant continues to op- 
erate fairly efficiently, while the totally-glued lubricant 
quickly loses its tribological properties. A similar behav- 
ior occurs for chains of length L = 3 and 7. 

The L-dependence of friction forces is illustrated in 
Fig. [9] for a configuration with exactly two molecular lay- 
ers {Qm = 2). In detail, the critical velocity Vc marks the 
crossover illustrated in Fig. [21 and the frictional forces js 
and /fc are evaluated immediately below and above Vc-, in 
the stick-slip motion and smooth sliding regimes respec- 
tively. Note that for the head-glued lubricant, the critical 
velocity depends only weakly on the molecular length L, 
while in the totally-glued case, Vc increases rapidly with 
L, so that for long molecules, L > 5, the critical velocity 
of the totally-glued chains is much higher than for the 
head-glued ones. Note also that the static friction /s for 
the stick-slip motion at ^ Vc is approximately twice 
as large as the kinetic friction fk in the smooth-sliding 
regime at Vg ^ Vc- 

We come now to discuss in some detail how the tribo- 
logical properties vary with the system parameters. Con- 
sider first the role of temperature. Figure [TUl reports the 
T-dependence of kinetic friction at a fixed driving veloc- 
ity Vs — 0.316, in the smooth sliding regime. Friction 
generally decreases when temperature grows, as is to be 
expected of the viscosity decrease of a fluid [l3|- Gen- 
eral trends for different lubricant molecules are similar: 
the head-glued lubricant provides the largest friction, the 
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FIG. 10: (Color online) Temperature dependence of the ki- 
netic friction force fk of the monoatomic lubricant {L = 1, 
N = 357, M2 = 100, diamonds) and of head-glued (full 
circles) and totally-glued (open triangles) chain molecules of 
length L = 7 (M = 51, iV = 357, M2 = 14), at fixed driving 
velocity Vs = 0.316. The inset shows the effective thermal 
expansion of the mean lubricant thickness. Error bars are 
smaller than the symbol size. Other simulation parameters 
same as in Fig. [3] 



monoatomic lubricant the lowest friction. It is interest- 
ing that the temperature dependence of kinetic friction 
fk for the totally-glued lubricant is significantly weaker 
than for the head-glued and monoatomic lubricants. The 
same type of dependence was observed for different driv- 
ing velocities. The monoatomic lubricant friction actu- 
ally increases with T, at low temperature, where the lu- 
bricant structure is crystalline. This is expected for an 
incommensurate lubricant/substrate interface [3]. The 
thermal expansion of the simulated lubricant is rather 
large, especially for the monoatomic lubricant, as shown 
in the inset of Fig. [TOl 

Next, Fig.fTTlshows a typical load dependence of the ki- 
netic friction. The friction force increases with load, but 
not quite linearly, i.e., the friction coefficient fj,k = fk/fi 
is not constant (unlike Amontons first law): fik decreases 
gradually the load /; increases. The reason is most likely 
associated with the increasing effective shearing temper- 
ature associated with a fixed shearing velocity under 
increasing lubricant compression, the lower viscosity in 
turn reducing friction. Note that friction of the head- 
glued lubricant is approximately twice that of the totally- 
glued system. 

Figure [T^ reports the dependence of the friction forces 
fs and fk on the driving velocity at a nonzero temper- 
ature T = 0.1. The friction forces generally rise with ve- 
locity, but very gently so. By comparison with Fig. [51 we 
note that the transition from stick-slip to smooth sliding 
at a nonzero temperature and for more than two layers 
occurs earlier than at T = and for exactly two layers. 
The transition shows a nonmonotonic dependence on the 
molecular length L, stick-slip persisting to the highest 
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FIG. 11: (Color online) Dependence of the kinetic friction 
force /fc and the friction coefHcient /x^ = jhl h (inset) on the 
load force for the totally-glued (dotted curve and open dia- 
monds) and head-glued (solid curve and diamonds) lubricant. 
Data refer to M = 70 chains of L = 5 atoms, other parameters 
same as in Fig. [S] Note the semilog scale. 
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FIG. 12: (Color online) Friction forces fs and fk as func- 
tions of the driving velocity for linear lubricant molecules of 
different lengths L — 1 [M = N = 360, solid diamonds), 
L = 3 (M = 120 and = 360, red circles), L = 5 (M = 70 
and A'' — 350, magenta up triangles), and L — 7 {M = 51 
and A'' = 357, blue down triangles), for either totally- (open 
symbols and dotted curves) or head-glued (solid symbols and 
curves) molecules. The calculation is carried out at tempera- 
ture T = 0.1, other parameters same as in Fig. [3] 



critical velocity for L ~ b. In the low-w^ regime the 
friction data tend to become noisy, announcing the well- 
known difficulty of MD to address the low shear rates of 
realistic lab experiments. For an extension of the present 
study to that reg ime, a transient-time correlation func- 
tion approach fl3| would be of substantial help. 

Figure [T3] demonstrates the dependence of the ki- 
netic friction on the strength of the interaction Vu be- 
tween the lubricant molecules. At small Vu, e.g., for 



FIG. 13: (Color online) Kinetic friction force fk as a function 
of the energy amplitude Vu of interaction between the lubri- 
cant atoms, for linear lubricant molecules of different lengths 
L = 1 {M ^ N = 360, solid diamonds), L = 3 (M = 120 
and N = 360, red circles), L = 5 (A/ = 70 and iV = 350, ma- 
genta up triangles), and L — 7 {M — 51 and A'^ = 357, blue 
down triangles), for either totally- (open symbols and dot- 
ted curves) or head-glued (solid symbols and curves) cases. 
Simulations carried out at temperature T = 0.1 and driving 
velocity Us = 0.1; other parameters same as in Fig. [3] 



Vu = 0.1, the surfaces are covered by monolayers of lubri- 
cant molecules, and the sliding interface lies somewhere 
near the middle of the lubricant film. As Vu increases, 
sliding becomes more and more viscous, thus friction in- 
creases reaching a maximum near Vu = ^ V^i =5. 
At this point, the sliding mechanism changes: the lubri- 
cant becomes rigid enough to move as a whole body, and 
frictional sliding is now shifted to the substrate/lubricant 
interfaces. If Vu increases further beyond this point, slid- 
ing becomes easier and easier because of increasing lu- 
bricant rigidity; this regime corresponds to the sliding of 
two hard incommensurate surfaces P, 0] ■ We found the 
same type of crossover for other numbers of lubricant 
molecules. 

Finally, Fig. [HI demonstrates the variation of the fric- 
tion coefficient with the substrate lattice constant Rs . In 
these simulations, we keep the size of the system con- 
stant at Lx ~ Q2 X 2/3 « 41.3, and vary the num- 
ber of substrate atoms from Ng = 41 to 83, so that 
Rs = Lx/Ns changes from « 0.5 to « 1. The total load 
force Fi = fiNg is also kept constant at the level Fi — Q2 
(so that we have fi = 1 for Ai"s = 62 as in the previous 
simulations). As expected, there appears to be no sys- 
tematic dependence of friction on the substrate/lubricant 
interface Rs/Ru for the head-glued system, where the 
molecules arrange themselves at frequently varying mu- 
tual distance, and no commensurability effect can arise. 

For the totally glued case, however, the friction coef- 
ficient shows a well pronounced minimum at Rs ^ 0.75, 
especially for a lower value of the lubricant-substrate in- 
teraction Vsi = 1. This is associated to a slow "cater- 
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FIG. 14; (Color online) Friction coefficient fi = ft / fi sls a, 
function of the substrate lattice constant Rs for AI — 70 linear 
lubricant molecules of length L = 5 for totally glued (blue 
open symbols) and head-glued (green solid symbols) chains. 
The sliding velocity is Vs = 0.316, the interaction radius R^i — 
Rs is varied, and in (b) all other parameters are as in Fig. |3] 
including Vai = 5 or 1/9, while in (a) the gluing interaction is 
Vs, = 1 or 1/9. 
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FIG. 15: (Color online) The central region of five successive 
snapshots of the motion of M = 70 totally-glued L = 5 lubri- 
cant molecules. The time interval between successive snap- 
shots (a) to (e) is At = 1. Snapshots b-e display only the 
atoms nearest to the bottom substrate, for better visibility of 
the caterpillar motion. To provide a clearer reference frame, 
one of the static atoms in the bottom-substrate is drawn as 
empty. The other parameters are as in Fig. I14f a') , except for 
Rsi ~ Rs = 2/3. 



pillar" motion of the glued molecules in the periodic 
substrate potential [l3|, caused by advancing of kinks 
(misfit dislocations) at the lubricant-substrate interfaces 
[Tsl [l^ . For a rigid linear molecule in a sinusoidal 
potential, virtually free sliding of the molecule should 
be observed for ratios Rs/Rnw\ — L/{2L — n), where 
n = 1,2, . . . ,int[L/2]. This gives, e.g., i?s/i?moi = 2/3, 
or Rs « 0.57 for L = 2 and 4, and Rs/Rmoi = 5/8, or 
Rs ~ 0.535 for L — 5. For a nonrigid molecule, these 
values are shifted upwards [l^ . As shown in Fig. [13 a 
caterpillar-type motion indeed takes place. However, as 
some molecules are attached with all L atoms but others 



stick only partially, e.g., with 2 — 4 atoms, while the glued 
layer should move as a whole (the molecules cannot over- 
take one another), caterpillar motion is disturbed, and 
this makes the minimum in Fig. [TJ] broad and shallow. 
The dependence of friction on the incommensurability 
between the substrate and molecular spacings was actu- 
ally observed experimentally [l8| . 



V. DISCUSSION AND CONCLUSION 

We presented simulations of a simple 2D model of 
molecular lubrication. Results indicate that layers of sim- 
ilar thickness of totally-glued lubricant molecules gener- 
ally produce a smaller friction than lubricant molecules 
that are head-glued. Totally-glued lubricants also stick 
better to the substrates, and boundary layers thereof may 
accordingly be harder to squeeze out under high load. 
On the other hand, whenever the applied load becomes 
strong enough to squeeze the lubricant out to less than 
two molecular layers, then the head-glued kind behaves 
much better in sustaining smooth sliding and low fric- 
tion even at very small residual coverage. We predict that 
when such extreme-load situations are likely to occur in a 
practical lubricated sliding devices, head-glued molecules 
could provide the best lubricant, unless the squeeze-out 
resistance is much worse than for totally glued lubricant. 
Moreover, when the base lubricant is totally- glued as a 
typical oil and, therefore, provides low friction proper- 
ties, then even a small concentration of head-glued ad- 
ditives such as diblock polymers may prevent the sub- 
strates from wearing at places (asperities) where the sur- 
faces come close to each other within a thickness of less 
than two layers. This is a new result, susceptible of ex- 
perimental test. 

Similarly, friction shows a clear increasing trend as a 
function of molecular length - the number L of monomer 
units in the molecules. This result is similar to the L de- 
pendency of the bulk viscosity, and our results compare 
well to those of Ref. [1]. This would make monoatomic 
liquids the best lubricants, if it was not for their much 
reduced resistance against squeeze-out, which eventu- 
ally makes polyatomic chains more effective lubricants 
at loads where boundary lubrication matters. 

The detailed search for the best compromise between 
the necessity of a squeeze-out resistance, calling for large 
L, and of low-friction characteristics, calling for small L, 
is a rather intricate task that is probably best undertaken 
on a case-by-case base by experimental trial and error 
procedures, or possibly by very extended simulations of 
sophisticated 3D realistic models. The simple 2D model 
at hand allows us to understand the main mechanisms 
at play in the soft lubricant and to predict the general 
trends. For example, local extra heating induced by the 
increased lubricant shear velocity at thinnest lubricant 
regions between substrate corrugations may well act as an 
important mechanisms producing a considerable overall 
friction reduction in macroscopic sliding. 
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